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Abstract. This paper presents a geometric-variational approach to contin- 
uous and d iscrete second-order field theories following the methodology of 
[MPS 1998 1 . Staying entirely in the Lagrangian framework and letting Y de- 



note the configuration fiber bundle, we show that both the multisymplectic 
structure on J 3 Y as well as the Noether theorem arise from the first varia- 
tion of the action function. We gener alize the m ultisymplectic form formula 
derived for first order field theories in [V1PS 199£], to the case of second-order 



field theories, and we apply our theory to the Camassa-Holm (CH) equation 
in both the continuous and discrete settings. Our discretization produces a 
multisymplectic-momentum integrator, a generalization of the Moser-Veselov 
rigid body algorithm to the setting of nonlinear PDEs with second order La- 
grangians. 
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1. Introduction 



This paper continues the development of the variational approach to multisym- 
plectic field theory introduced in Marsden, Patrick, and Shkollcr [ MPS 1998| . In 



that paper, only first-order field theories were considered. Herein, we shall focus 
on second-order field theories, i.e., those field theories governed by Lagrangians 
that depend on the spacetime location, the field, and its first and second partial 
derivatives. 

Multisymplectic geometry and its applications to covariant field theory and non- 
linear partial differential equations (PDE) has a rich and interesting history that 



we shall not discuss in this paper; rather, we refer the reader to |Ga 1974, G 1991 



GIM 1999, MPS 1998 1 and the references therein. The covariant multisymplectic 
approach is the field-theoretic generalization of the symplectic approach to classi- 
cal mechanics. The configuration manifold Q of classical Lagrangian mechanics is 
replaced by a fiber bundle Y — > X over the n+l dimensional spacetime manifold 
X, whose sections are the physical fields of interest; the Lagrangian phase space is 
TQ in Lagrangian mechanics, whereas for kth-order field theories, the role of phase 
space is played by the kth-jet bundle of Y, J k Y, thus reflecting the additional 
dependence of the fields on spatial variables. 

For a given smooth Lagrangian L : TQ — > R, there is a distinguished sym- 
plectic 2-form Wi on TQ, whose Hamiltonian vector field is the solution of the 
Euler-Lagrange equations of Lagrangian mechanics. Lagrangian field theories, on 
the other hand, governed by covariant Lagrangians C : J k Y — > A n+1 (X), can be 
completely described by the multisymplectic n+2-form on J 2k ~ 1 Y, the field- 
theoretic analogue of the symplectic 2-form of classical mechanics. In the case 
that X is one dimensional, tie reduces to the usual time-dependent 2-form of clas- 



sical nonautonomous mechanics (see [ MS 1999 1) 



Traditionally, the symplectic 2-form lol as well as the multisymplectic n+2-form 
flc are constructed on the Lagrangian side, using the pull-back by the Legendre 
transform of canonical differential forms on the dual or Hamiltonian side. Recently, 
however, Marsden, Patrick, and Shkoller [ MPS 1998| have shown that for first order 



field theories wherein C : J Y Y — > A n+1 (X), flc = dQc arises as the boundary term 
in the first variation of the action J x Coj 1 ^ for smooth mappings <f> '■ X — > Y. This 
method is advantageous to the traditional approach in that 

(i) a complete geometric theory can be derived while staying entirely on the 

Lagrangian side, and 
(ii) multisymplectic structure can be obtained in non-standard settings such as 
discrete field theory. 

to the 



The purpose of this paper, is to generalize the results of [MPS 1998 



case that C : J 2 Y — > A n+1 (X). In Section 2, we prove in Theorem 2.1, that a 
unique multisymplectic n+2 form arises as the boundary term of the first variation 



of the action function. We then prove in Theorem 2.2, the multisymplectic form 
formula for second-order field theories, a covariant generalization of the fact that 
in conservative mechanics, the flow preserves the symplectic structure. We then 
obtain the covariant Noether theorem for second-order field theories, by taking the 
first variation of the action function, restricted to the space of solutions of the 
covariant Euler-Lagrange equations. 

In Section 3, we use our abstract geometric theory on the Camassa-Holm (CH) 
equation, a model of shallow water waves that simultaneously exhibits solitary 
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wave interaction and wave-breaking. We show that the multisymplectic form for- 
mula produces a new conservation law ideally suited to study wave instability, and 
connect our intrinsic theory with Bridges' theory of multisymplectic structures (see 



B 1997| and flMPS 1998| ) 



Section 4 is devoted to the discretization of second-order field theories. We 
are able to use our general theory to produce numerical algorithms for nonlinear 
PDE that are governed by second-order Lagrangians which naturally respect a 
discrete multisymplectic form formula and a discrete Noether theorem. Again, we 
demonstrate this methodology on the CH equation. 

2. Variational Principles for second order Classical Field Theory 

2.1. Multisymplectic geometry. In this section, we review some aspects of mul 
tisymplectic geometry following |G 1991[ , JGIM 1999[ , ]MS 1999| , and ]MPS 1998 1 



Let X be an orientable (n+l)-dimensional manifold (which in applications is 
usually spacetime) and let ttxy ■ Y — > X be a fiber bundle over X. Sections 
<j> : X — ► Y of this covariant configuration bundle will be the physical fields. The 
space of sections of ttxy will be denoted C°°(itxy) or by C°°(Y). The vertical 
bundle VY is the subbundle keiTwxY of TY, where Tttxy denotes the tangent 
map of the ttxy- 

If X has local coordinates x^, [i = 1, 2, . . . , n, 0, adapted coordinates on Y are 
y A , A = 1, . . . , N, along the fibers Y x := tt xy (x), where x & X and N is the fiber 
dimension of Y . 

J k Y denotes the k th jet bundle of Y, and this bundle may be defined inductively 
by J x (- ■ ■ (^Y)). Recall that the first jet bundle J Y Y is the affine bundle over 
Y whose fiber over y G Y x consists of those linear mappings 7 : T X X — > T y Y 
satisfying 

Tttxy 07 = Identity on T X X . 

Coordinates (x^, y A ) on ttxy induce coordinates y A on the fibers of J X Y. Given 
4> G C°°(Y), its tangent map at x G X, denoted T x cf>, is an element of J 1 !^^). 
Therefore, the map x 1 — ► T x (j) defines a section of J X Y regarded as a bundle 
over X. This section is denoted j 1 (4>) and is called the first jet of </>, or the first 
prolongation of (p. In coordinates, j 1 ((f>) is given by 

x^^(x^,cf> A (x^,d u( f> A (x' i )), 

where d v = d/dx u . A section of the bundle J X Y — ► X which is the first prolon- 
gation of the section of Y — ► X is said to be holonomic. 

The first jet bundle J X Y is the appropriate configuration bundle for first-order 
field theories, i.e., field theories governed by Lagrangians which only depend on the 
spacetime position, the field, and the first partial derivatives of the field. Herein, we 
shall focus on second-order field theories that are governed by Lagrangians which 
additionally depend on the second partial derivatives of the fields; thus, in second- 
order field theories, the Lagrangian is defined on J 2 Y = J 1 (J 1 F). Let us be more 
specific. 

Definition 2.1. The second jet bundle is the affine bundle over J X Y whose fiber 
at 7 G J x Y y consists of linear mappings s : T X X — ► T 7 J X Y satisfying 

Tttx,,j 1 y s — Identity on T X X . 
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One can define the second jet prolongation of a section <j> : X — ► Y, j 2 (4>), as 
that is a map x i — ► T x j x ((f)), where j 1 (<j>) is regarded as a section of 
over X . This map defines a section of J 2 Y regarded as a bundle over X with 
j 2 ((f>){x) being a linear map from T X X into Tji^\r x \ J X Y. In coordinates, j 2 (</>) is 
given by 

x» — » (x^^^),^^^),^^^^)). 

We shall also use the notation <^ lM2 = d^d^^ for second partial derivatives. A 
section p of J 2 y — ► X is said to be 2-holonomic if p — j 2 (^Y,.i 2 Y ° P)- Continuing 
inductively, one defines the k th jet prolongation of <f>, j k ((f>), as • ■ (j 1 (<)>)))■ 

Consider a second-order Lagrangian density defined as a fiber-preserving map 
C : J 2 Y — > X n+1 {X), where X n+1 (X) is the bundle of (n+l)-forms on X. In 
coordinates, we write 

£(s)=L(x», y A ,yl,y$ lfi2 )u>, 

where ui — dx 1 A ... A dx n A dx°. 

For any fcth-order Lagrangian field theory, the fundamental geometric structure 



is the Cartan form 0£; this is an (n + l)-form defined on J 2fc ~ 1 Y" (see [G 1991]). 
For second-order field theories, the Cartan form is defined on J 3 Y, the covariant 
analogue of the phase space in mechanics. The Euler-Lagrange equations may be 
written intrinsically as 

(j 3 (f>y(V J dS c ) = V V g T(J 3 Y), (2.1) 

where J denotes the interior product. Traditionally, the Cartan form is defined 
using the pull-back by the covariant Legendre transform of the canonical multisym- 



plectic n+l-form on the affine dual of J 2k ~ x Y (sec [ p 1991| , |GIM 199S| , |MS 1999Q 
In local coordinates, the Cartan form on J 3 Y is given by 

( dL _ ( dL \\ , A dL A 

dL a „ / dL \ A dL 



L -^~ +D '{M.) y "-^ ,: ") u - (2 - 2) 

where uj v = d v J uj and uj^ — d v J J w, etc. For a fc th -order function / e 
C°°( J k Y, K), the formal partial derivative of / in the direction x^, denoted by 
D^f, is defined by (j k+1 (f>)* (D^f) = o j k <j)) for all <j) 6 C°°(Y), and is a 

smooth function on ,] k+1 Y . In jet charts, 



D„ f ( -w + ^ + - + ^-j/ / - 1 ,.., I*.* 



In the next section, we shall prove that the Cartan form arises as the boundary 
term in the Lagrangian variational principle. 



2.2. The variational route to the multisymplectic form. In this subsection 
we show that a multisymplectic structure is obtained by taking the derivative of an 
action functional, and use this structure to prove the multisymplectic counterpart 
of the fact that in conservative mechanics, the flow of a mechanical system consists 
of symplectic maps. 
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Let U be a smooth manifold with (piecewise) smooth closed boundary. Define 
the set of smooth maps 

C°° = {</> : U — ► Y ttxy °4>-U — > X is an embedding} . 

For each <j) G C°° set <fix '■= t^xy 4> an d Ux 4>x(U) so that <\>x ■ U — ► Ux is 
a diffeomorphism. Let C denote the closure of C°° in some Hilbcrt or Banach space 
norm. The choice of topology is not crucial in this paper, and one may assume all 
fields are smooth. The tangent space to the manifold C at a point <p G C is given by 

{V G C°°(X, TY) | tty^ty o V = and Vx := T7r xy ol/o^isa vector field on X} . 

Consider G, the Lie group of 7rxy-bundle automorphisms i]y ■ Y — ► Y covering 
diffeomorphisms r)x ■ X > X. 

Definition 2.2. The group action $ : G x C — ► C is given by 

${nY,4>) =Vy of 

Note that (jiy°4>)x — Vx°<f>x, and if<f>o<j> x G C°°(ttu x ^y), then (r]Y (t>) <f>x or 1x G 

The fundamental problem of the classical calculus of variations is to cxtrcmizc 
the action functional over the space of sections of Y — ► X. 



Definition 2.3. The action functional S : C 



is given by 



S(4>) 



JUx 



(2.4) 



Definition 2.4. (p G C is said to be an extremum of S if 



d_ 

dX 



S(*(f£,0))=O 



for all smooth paths Ahi)| in G, where for each X, n Y covers n x . 

One may associate to each 4> x G C, the section of Y given by f^o^o^'jofr)^)" 1 ; 
namely rj Y ° (4> ° ^jf 1 ) ° (Vx)^ 1 ma P s ^x : ~ ° 4>x(U) into <p x (U). 

If we choose the curve ^ A such that <j>° = <f> and (d/dA)lA=n$(?7v-, 4>) = V, then 
we have that V 
following: 



(d/dA)| A=0 <^ andT/ x - £| x=0 »&. 



This will be used in the 



d<S -V 



d_ 

d 

dA 
/ - 



d 



£(i 2 (^ A o(^)- 1 ))+ f 



d 



dA 



(^)*£( J 2 (0o^- 1 )) 



E/x 



iVxA/OWx)) , (2.5) 



where * stands for the pull-back, and £ denotes the Lie derivative. 

Now, let VY C TY be the vertical subbundle; this is the bundle over Y whose 
fibers are given by 



V y Y = {ve TyY I Titxy ■ v = 0} . 
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For each 7 £ J Y Y y there exists a natural splitting T y Y = image 7 © V y Y. For 
example, for a vector V £ T^C, let 7 = fty o fa 1 ), V h := y(V x ), and V v := 
Fo^ 1 - V'\ Then 

Titxy o V h = Titxy ° t(Vx) = idrx(Vx) = Vx ■ 

On the other hand, by definition, Vx — Tttxy °V o <f> x . Therefore, Tttxy -V v — 
which confirms that any vector V £ T^C may be decomposed into its horizontal 
component 

V h = T(4>o#£)-Vx ) (2.6) 

and its vertical component 

V v = V o fa 1 - V h . (2.7) 

Remark 2.1. Notice that V[x) £ T^ X )Y for all x £ U, while V h and V v are 
vector fields on Ux = 4>x{U). 

Next, we define prolongations of automorphisms rjy of Y and of elements V £ 
T^C. 

Definition 2.5. Given an automorphism rjy of Y — » X , its first prolongation 
^{Vy) ■ J X Y — ► J l Y is defined via 

3 X {t]y){i) = Tr lY o 7 o Tri x x . 

If 7 : T X X — > T y Y, then j l {riY){l) ■ T VX / X \X — ► T VY / y -\Y, with local coordinate 
expression 

To define the first prolongation of a vector V £ T^C, denoted J 1 (l^), let r)y be a 
flow of a vector field v on Y with v o <fi = V. 

Definition 2.6. The first prolongation ^{V) of V is a vector field on J X Y given 
by 



jHvy) 



If in a coordinate chart V = (V^, V A ); identifying V with V o (fix 1 , we see that 
( p.8| ) yields the following local expression for j 1 (V^)(7): 

Using induction, one can define the fc th prolongation of an automorphism Tjy and 
the k th prolongation of a vector V £ T^C for all k > 1, and these will be denoted 
by j k { 7 ]Y) and j k (V), respectively. 

Definition 2.7. For a k th - order function f £ C°°(J k Y,M.), the variational deriv- 
ative of f is the function on J 2k Y given by 



1L-^ (_i YD ...n ( d f 
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In particular, for a second order function / G C°°(J 2 Y,K), the variational de- 
rivative of / is the function on J 4 Y given by 

5f - df dJX\+d„d( 01 



Sy A dy A v \dy A ) v M \dy A 
Definition 2.8. Let C 4 = {j 4 (<j) o <jy^)\<j) G C}. 

Theorem 2.1. Given a smooth Lagrangian density C : J 2 Y — > A n+ (X), there 
exist a unique ^ G K n+2 {J 4 Y) given by 

5L 

Sy J 

a unique map V EL C G C°°(C 4 , T*C <g> A n+1 (X)) given by 

V EL £(<f>).V = j\ ( t>o<t>x 1 r(J^iv(dy A Au>)} , (2.10) 

and a unique differential form 0£ G A n+1 (J 3 Y) given by 

( dL _ ( dL \\ , A dL A 

dL A „ { dL \ A dL 



L -w^ +D 'WJ y: ~^ ,: ')" (2 - n) 

swc/i that j 3 ((f> 4>x l )*®C = £ ° J 2 (<^ ° '/'x 1 ) / or an 2/ ^ ^ and £/ie variation of 
the action functional S is expressed by the following formula: for any V G T^C and 
any open subset Ux of X such that Ux f] dX = 0, 

dS^-V= f V EL C{<p)-V+( f^o^nf^JQc}. (2.12) 

JUx JdUx 

Furthermore, 

V EL C(cj ) )-V = j 3 (cj ) o^~ 1 )*[f(V)jnc} inU x , (2.13) 
where Q,£ = d0£ is the multisymplectic form on J 3 Y . The variational principle 



(2.1i) yields the Euler-Lagrange equations l[2.y on the interior of the domain, 
which in coordinates are given by 

° *?» " {wi^ °^») + ^ (H^(* • «'») = ' 

(2.14) 

while the form 0£ naturally arises in the boundary term and matches the definition 
of the Cartan form given in (2.i). 

Proof. The proof proceeds in three steps. We begin by computing the first variation 
using (2.5). Then we show that the boundary term yields the Cartan form. Lastly, 
we verify the statements related to the interior integral. 

Choose Ux — 4>x{U) small enough so that it is contained in a coordinate chart. 
If in these coordinates V — (V^, V A ), then along (j}0(j) x 1 , the coordinate expressions 
for Vx, V , V v are written as 

Vx = v—, v h = v»-?- + v» di<l>0<l> x 1)A 9 



dx^ ' dx^ dx^ dy A ' 
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and V" — (V v ) — — r := [ V — 

QyA 



,d{4>o(t> 



-1\A 



X 



dx^ 



W ' (2-15) 



Using the Cartan formula we first compute the second term on the right-hand- 
side of Q: 



£ Vx C(j 2 (cj ) o ( j >x 1 )) = / £v x (Lu) 



U x 



di Vx (Luj) + i Vx d{Luj) 



DUx 



Li Vx uj = 



L V ti oj e . 



dU x 



(2.16) 



Using (2.7), and the local expression for the vertical vector field V v , we have 
that 



Ux 



d_ 

dX 



£(iV°(^r x )) = 



dL 



U x 



dL 



dL 



(2.17) 



In the following, we shall use D v f for the formal partial derivative of a function 
/ (see (2.3)), and will denote gfr(/ j 2 (<fr ° ^x 1 ))- Integrating ( 2.17 ) by parts, 
we obtain that 





' dL 


- — { 


/ 

'Ux 


dy A 


dx v \ 



dL 



{V V ) A LU 



CO ■ 



Ux 



dL 



, \9y A 



(vXu " 



( d f dL 



Ux \ dxt " \ d V A 



(V 



v\A 



Using the fact f >u u> — d(fu> v ), applying the Stoke's formula J v da = J gu a, and 
combining the last calculation with ( 2.16] ), we obtain 



dS A ■ V = 



Ux 



' dL 


d 


(8L\ 


d 2 


dy A 


~ dx" 


\dy A ) 


dxvdx^ 



( dL 



(V v ) a uj 



(2.18) 



dL d f dL 



dT 

(O V + #x(v^)X + L ye ^ ■ 



'dux \ d v A dxt * \ d y A p,, 

Definition 2.9. A form a on J k Y is contact, if (j k (f>)*a — for all 4> e C°°(Y). 



Lemma 2.1. For a smooth Lagrangian density C : J Y — > A n+ (X) there exist a 



unique diffe renti al form 0£ £ A n+1 (J 3 y) defined by ( 2.11 ) such that the boundary 
integral in (2.18) is equal to 



j 3 (<f>°<t>xT\3 3 (y)JQc]. 

dUx 

Furthermore, 0£ can be written as a sum of Luj and a linear combination of a 
system of contact forms on J 2 Y with coefficients being functions on J 3 Y. 
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Proof of Lemma Let W = (W^, W A , W A , W A U , W A v8 ) be an arbitrary vector 
field on J 3 Y, and let ip :— j 3 {<po(j)^}) 7 a map from Ux to J 3 Y. Then one computes 

\w{dy A A lo v ) = W a uj v - W 6 dy A A w v8 , 

i w (dy A A Up) = Wfun - W e dy A A w Mfl , 

8(6 o 6~ l ) A 

<p*i w (dy A A u)y) = W A u v - W e KV J x ' dx" A lj vB ■ 

ox^ 

Using the formula 

[ if/z^,0 
dx" A uj v e = I u u if fi = 9 (2.19) 
[ -oj0 if n = v, 

one finds that 

<P*iw(dy A A «„) = W A u v - W e w J* ' uj v + q * W % • 
Similarly, 

* , a a Rd 2 (d>o6' 1 ) A d 2 Uo6- 1 ) A fl 



Thus if we let W — j 3 (V), use (|2.9|) , and recall the local expression ( 2.15| ) for 
(V V ) A , we obtain that 

<p*i j3(v) (dy A A w „) = (V*)V + d( *° d ** )A V 6 ue , 

^i i3(y) (^ A w „) = (Vtf ^ + . 
Next, observe that V wj = iyw. Also, 

3 J* = A* ° *x)V. and = f{<P ° • 

These observations together with the previous identities imply the following impor- 
tant formulas 

j S (<f>°<t>xy\j 3 (V)j(dy A A^-y A ^)} = (V")>„, , . 

J 3 (0 o ^)*[7 3 (V) J (dy A Acv v — y A u)\ = (V«) a uj v . l ^ Uj 

Subst ituting these formulas into the boundary integral of the variational principle 
( 2. IS ) , we obtain that 

dL d ( dL 



du x \ d V A dxt " \ d U, 



ftt o &r\j 3 (V)j \(^- D J^)) (dy A - 



9y^ 



i\dy A Dtl i. dy»n ) ) dy AuJl/ + dy A ^ dVv 



T dL A ( dL \ A dL A \ 
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JdUx 

This proves the existence of a unique differential form 0£ and demonstrates how 
this form naturally arises in the boundary integral of the variational principle. 
Integration by parts yields the boundary integral with terms that involve partial 
derivatives of (V V ) A of all orders up to k — 1 (in our case k — 2). Equation (2.20) 
show that each partial derivative of (V V ) A has an associated (n + l)-form on J 2 Y, 
and substitution of these forms yields a unique differential (n + l)-form as desired. 

Since L and its partial derivatives are functions on J 2 Y, then by (2.3), ( gpi ) 

is a function on J 3 Y, and therefore 0£ is a (n + l)-form on J 3 Y. 
It is easy to show that 

j k ^o^ x 1 )*(dy A Au; ly -y A Lo)=0 

j k (<j ) ocj) x 1 y(dy A Au;,-y A ^)=Q 

for all integers k > 2 and for all <j> G C. Therefore, dy A f\u) v —y A u) and dy A Aujfj,—y A u> 
are contact forms on J 2 Y. Hence the last statement of the lemma follows. □ 
A simple computation then verifies that 0£ is the Cartan form so that 

/(0o0- 1 )*e £ = £ oj 2 (^otf). 

Next, consider the interior integral of the variational principle ( 2.1§| ). Since 
j k ((j) o (f)x 1 )*ijk/ v \(dy A Aw) = (V v ) a uj for all integers k > 1, we obtain that 





' dL 


d 


(8L\ 


d 2 


/ 

lUx 


dy A 


dx v 


\dy A J 


dx v dx^ 



( dL 



l'x 



3 i {4> 4> x 1 )*' l o i {V) 



<°4>x l ) 



U x 



dL 
dy A 
5L 



dL 



rv\A, 



dy A 



dL 
dy A n 



dy A A lj 



-dy A Aw 



(2.21) 



where Jpr is the variational derivative of L in the direction y A (see Definition 2.7). 
Since L is a function of second order by hypothesis, then its variational derivative 
is a function on J 4 F. Therefore the form 'J = j^dy A A a; is an (n + 2)-form on 



SL 



iv{dy A A to) 



J 4 Y. Moreover, the integrand in (2.21) written as j 4 (<j) ° fix 

defines a unique smooth section V EL C G C oc {C i 1 T*C® A n+1 {X)) as desired in the 
statement of the theorem. 

Now we shall prove the following 

Lemma 2.2. The forms flc = d&c and VP = jp?dy A A u) satisfy the following 
relationship: 



/(0 O Jf 1 )*i J 4 (V) * = /(0 O Jf 1 )*i J 3 (y) 



CI i 



(2.22) 



for all 4> G C and all vectors V G TC. 

Furthermore, a necessary condition for <fi G C to be an extremum of the action 
functional S is that 

f(^o^ x 1 )*i w n c = (2.23) 
for all vector fields W on J 3 Y , which is equivalent to 



J 



'x 



= 



(2.24) 
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for all vector fields V on J A Y . 



Proof of Lemma 2.L The proof will involve some lengthy computations that we 
partially present below. To compute flc-, let us write 0£ as 



Or = 



dy A D " 



dL 
9y± 



{dy A Au v - y A u) 



dL 



Then, for W G TJ 3 Y, we obtain 



4 {dy A A^- y^uj) + Luj. 



i w Q c = W 

- d 
+ W 

- d 



' dL 

dL 



dL 



- D 
D 



dL 



M \dy A , 
dL 



(dy A Aw„ - y A u) 

A {W A u v - W 9 dy A A u v9 - y A W e ujg) 



dy»n 



(dy A A - y^ui) 



D, 



( dL 



\9y± 



FIT 

(W A to - W e dy A A we) + ^a( wA u ~ W e dy A A uj e ) 



The last step is to pull back iw^c by v 5 = j 3 ( < /' ° ^x 1 )' ^ n is eliminates the terms 
with the contact forms. In addition, using the fact that the pull back commutes 
with the exterior derivative, and applying formulas such as (2.19), we obtain that 



<p*i w n c = W J 



dL 

dy A 



d f dL 



dL 



dy A dx» \dy A ^ 



A u) v 



°4>x X ) A , A §_ ( dL \ d 2 ^- 1 

Kdx" we j dx^ \dy~j^J dx"dx" 

dL 



dx v 



-UJ0 



9(<f>°4> jl 

~ r dx" 
i " dxv 

v \dx» \dy A ^ 



dx" f 

1\A 



-UJ0 



dL 9{4>°<t> > 



Oy- 



dx^ 



dyt 



Aw,, 



Some cancellation and further rearrangement yields 



ip*i w £l c = W A 
- W 



dL 



d f dL 



dy A dx v \dy A 

gd^ocj)- 1 ^ ( dL 



d 2 



( dL 



dx 9 



dxvdx^ ydy^ 
d ( dL\ d 2 



( dL 



dy A dx" \dy A J dx v dx^ \dy A 



Letting W = j 3 (V), we have that 



dL 



<p'i jHV )toc = (v v ) 



d f dL 



d 2 



dL 

dy A dx" \dy A J ' dx^dx^ \dy A 



where the right hand side equals j i ((j)0(j) x 1 )*iji^ v ^ by ( |2.21 ). Hence, the relation 
( p^ ) is proved. 
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A necessary condition for <f> £ C to be an extremum of the action functional 
S is that the interior integral in (2. IS) vanish for all vectors V £ TC. From the 
calculation above, one may readily see that it is equivalent to the condition (2.23). 

Now if we let V be a vector field on J 4 Y, then 

IyV|/ = iv {j^ dyAhUJ 

= j^j(iv(dy A ) A uj - dy A A i v w) 
= ^ x {V A u-V 6 dy A Au e ) . 

Hence, 



dL 



d 
~dxP 



,dy A 

Thus, the condition 



dL 

dy A 



d 2 



( dL 



dx v dx^ \dy A 



V — V 



e ° 0x ) 



l\A 



dx e 



/(fio^riy* = 



for all vector fields V on J Y is equivalent to the condition (2.23). This completes 
the proof of the lemma. □ 

Lemma 2.2 contains two equivalent conditions for cf> £ C to be extremal. Both 
conditions yield the same coordinate expression of the Euler-Lagrange equations 
given by 



dL 



A (i 2 0Wx)) 



d ( dL 2 x 



dy AW ^ " irA " dx v \dy AyJ vr " ^ x 
which is the final statement of the theorem 



d 2 



dx u dx> 1 \dy A 



dL 



4 (/W> °<tf)) 



o , 



□ 



Remark 2.2. As one may see the proof we have presented can be generalized to 
Lagrangian densities on J k Y . One has to modify the labeling of variables to reflect 
the general case. For example, 

(^)l... w Wfl = V*[7 8 (V)J(d<... w A^-Vi...^)} , 

where < I < (fc — 1). Then the Cartan form shall arise in the boundary integral 
as a linear combination of the forms above. 

We shall call critical points 4> of S, solutions of the Euler-Lagrange equations. 

Definition 2.10. We let 

V = {<j)EC | f((t>o<t> x 1 )*i w Q. c =Q for all vector fields Won J 3 Y} (2.25) 

denote the space of solutions of the Euler-Lagrange equations. 

We are now ready to prove the multisymplectic form formula, a covariant gen- 
eralization of the symplectic flow theorem to second-order field theories. R 



1 For first-order field theories, this is Theorem 4.1 in MPS 199£] 
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2.3. The multisymplectic form formula. If </> A is a smooth curve of solutions of 
the Euler-Lagrange equations in V (when such solutions exist), then differentiating 
with respect to A at A = will give a tangent vector V to the curve at <j> = <fp. By 
differentiating -A\. n j 



-)~ 1 )*{W J fl/r] = 0, we obtain 

j 3 (4>°4>x 1 y £ 3 hv)[w _in £ ] = o 

for all vector fields W on J 3 Y. Therefore if V is a submanifold of C, then for 
any <j> E V we may identify T$P with the set of vectors V that satisfy the above 
condition. However, we do not require V to be a submanifold. 

Definition 2.11. For any 4> e V , 

T = {V E T C | j 3 (</)o ( j)- 1 )*£ j 3 ( y ) [W J n £ ] = /or a// vector fields V on J 3 Y} 

(2.26) 

defines a set of solutions of the first variation equations of the Euler-Lagrange 
equations. 

Theorem 2.2. (Multisymplectic form formula) If<f>EV, then for all V and 
W m T, 



0. 



(2.27) 



dU x 



Proof. We follow Theorem 4.1 in |MPS 1998| ]. Define the 1-forms ai and a 2 on C 

by 



ai(<t>)-V 



J 

JUx 



and 



Oa(0) • V := 
so that by (|2.12|) and (|2.13|) , 



J 



OUx 



dSj, ■ V = ai{4>) ■ V + a 2 {<j>) ■ V for all VeT^C 



(2.28) 



Furthermore, 

d 2 S{(j)){V, W) = dai((j>)(V, W) + da u (<f>)(V, W) for all V, W E T^C . 

Since d 2 S = 0, we have that 

dai(0)(V, + da 2 (</>)(^, W) = for all V,W E T^C . (2.29) 

Given vectors V,W E T^C we may extend them to vector fields V, W on C by fixing 
vector fields u, u; on Y such that V = v o and = u> o 0, and letting V(p) — v o p 
and W(p) = w o p. If 77^- covering 77^- is the flow of u, then ^( rjy, p ) is the flow of 
V. Notice that V(<f>) = V and W((f>) — W, hence the equation ( 2.29D becomes 

rfai(V,W)(0) + da 2 (V,W)((/)) = . 

Recall that for any 1-form a on C and vector fields V, W on C, 

da(V, W) = V[a(W)] - W[a(V)] - a([V, W}) . (2.30) 

Also recall that for a ve ctor field V on C and a function / on C, V[/] = df ■ V. We 
now use the latter and ( 2.30| ) on a 2 . We have that 

da 2 (V, W)(<f>) = V[aa(W)]($ - W[a 2 (V)](<£) - a 2 ([V, W])(<£) 
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= [tf(a a (W))-V](0) - Ka 2 (V)) • -a 2 (0) ■ [V,W] 

= d(a 2 (W))((f>) ■ V - d(a 2 (V))(^) • VF - a 2 {4>) ■ [V, W] . (2.31) 

Similarly, 

da 1 ((f>){V,W)^d(ai(W))((f>)-V-d(a 1 {V)){<l))-W~a 1 {<l))-[V,W} . (2.32) 
Let (f> € V and <fi x = r\y o be a curve in C through such that 

d 



V = 



dX 



and V e T . 



Now we restrict V, W to T . We shall give a detailed computation of the first term 
on the right hand side of (2.31). We have that 



d(a 2 (W))(<P)-V = ± 

d_ 

~ dX 
d 

~ dX 

- A 

~ dX 



(a 2 (W))(cf> x ) 
a 2 {4> X ) ■ (wocj) X ) 

jV°(^)~W(w A )-ie* 

v=o J dUx 



dU, 



f^orxT£p { v)(f{W)AQc) 
f^o^ x 1 )*d[j\V)Af{W)AQ L 



+ / j"^o^y[^(v)jd(f(w)jec)] , 

JdUx 

where the last equality was obtained using Cartan's formula. We have also used 
the fact that W x — w o A and W = W o <fi have the same fc-th prolongation. 
Furthermore, using Stoke's theorem, noting that ddUx is empty, and applying 
Cartan's formula once again to d(j (W) J &£.), we obtain that 



d(a 2 (W))(<f>)-V = / j 6 {<i>o<l)-£T\j 6 {V)J£ 3 z { w ) Qc] 



dUx 



dUx 



(2.33) 



Similarly, 



d{a 2 {V)){<t>)-W 



dUx 



dUx 



f^o^ x 1 r[f(w)jf(v)jn c ] 



(2.34) 



Now, j 3 ([V, W]) = \j 3 (V),f(W)]; hence, 
a 2 (4>)-[V,W] -- 



dUx 



3 3 ^o^y{[f{V),f{W)]A@ c ) 
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Recall that for a differential form a on a manifold M and for vector fields X, Y on 
M, 



Therefore, 
a 2 (<t>) ■ W,W] 



dUx 



dUx 



f (^o^TU 3 (v)j£p (w) ec-.f(v)jf(w)jn c -j 3 (w)j £ P {v) q £ ], 

(2.35) 



where we have again used Stoke's theorem and Cartan's formula twice. Substituting 
( |jf ), (|2.34| ), and (|2.35D into ( |2.3lD , we obtain that 



da 2 mV,W)= / j s (cf ) o^- 1 )*\j 3 (W)Jj 3 (V)Jfl 



dUx 



(2.36) 



We now compute ( 2.32). Similar computations as above yield 



d( ai (W))(4>)-V 



Ux 



j 3 (4>o4> x 1 )*£ j 3 {v) (f(W)jQ j: ) 



which vanishes for all G V and V G T. Similarly, d(ai(V))((f) ■ W = for all 
4> G T 3 and W E J 7 . Finally, ai(^) = for all ^ e P. Therefore, the equation 
( 2.32) vanishes for all <f> G 7-* and V,W € J 7 . Using the latter and ( 2.36| ), equation 
( ^.29| ) becomes 



dUx 

for all <b G V and all V.W £ J- , as desired. 



□ 



2.4. The Noether Theorem. Suppose that S is invariant under the action &(g, <ft) 
of a Lie group G on C. This implies that for each g G G, $>(g,<j>) G T 5 whenever 
G P. We restrict the action to elements of P. For each element £ of the Lie 
algebra g of G, let £c be the corresponding infinitesimal generator on C restricted 
to elements of V. By the invariance of S, 

S($(exp(t£), 4>)) = S{c/)) for all t . 

Differentiating with respect to t at t = 0, and using the fundamental property of 
the Cartan form that C o j 2 ((j> o (fix 1 ) = J 3 (0 ° ^x 1 )*®^' we ^ no - * na ^ 

(7x 



Then by Theorem 2.1 and the invariance of S we have that 



j J (<Wx TLf (&(</>)) J fi/ 



(2.37) 
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Definition 2.12. Let J G Hom(Q, T*C ® A"( J 3 Y)) satisfy 

j 3 (e c (^))jr! £ = d[j(O(0)] 

for all £ G g and (f> £ C. Then the map JJ : C — ► g* defined by 
zs i/ie covariant momentum map of the action. 



(2.38) 
(2.39) 



With this definition, ( 2.37 ) becomes J Ux d[j 3 ((j) o f^ 1 )* (J(</>), 0] = 0) an d since 
this holds for any Ux C X, the integrand must also vanish; thus, 

4f(<Wx)*WU>]=0. (2.40) 
On the other hand, by Stoke's theorem we may also conclude that 

j^o^rw),o=o. 



(2.41) 



Last two statements are equivalent, and we refer to them as the covariant Noether 
Theorem. 

3. A MULTISYMPLECTIC APPROACH TO THE CAMASSA-HOLM EQUATION 

3.1. The Camassa-Holm equation. The completely integrable bi-Hamiltonian 
Camassa-Holm (CH) equation^ 

Ut - Uyy t = SltUy + lUyUyy + UUyyy (3.1) 

is a model for breaking shallow water waves that admits peaked solitary traveling 
waves as solutions (see [ CH 1993 ], [ CHH 1994 1 ) . Such solutions, termed peakons, 
develop from any initial data with sufficiently negative slope, and because of the 
discontinuities in the first derivative, these solutions are difficult to numerically 



simulate, particularly in the case of a peakon-antipeakon collision (see [CHH 1994 



The multisymplectic framework for the CH equation is intended to provide 
a foundation for numerical discretization schemes that preserve the Hamiltonian 
structure of this model, even at the discrete level. After developing the multi- 



symplectic framework for (3.1), we shall follow [MPS 1998| and develop the entire 
discrete multisymplectic approach to second-order field theories, concentrating on 
the discrete CH equation as our model problem. Although we shall only produce 
the simplest multisymplectic-momentum conserving algorithm for this equation, our 
construction is completely general and will allow for the creation of fcth-order ac- 
curate schemes for ar bitr arily large k. 

The CH equation (3.1) is usually expressed in terms of the Eulerian, or spatial 
velocity field u(t, y), and is the Euler-Poincare equation for the reduced Lagrangian 

l(u) = \ f (u 2 + ul)dy. (3.2) 



Alternatively, one may express (3.1) in terms of the Lagrangian variable rj(t,x) 
arising from the solution of 

d 

— r){x,t) = u(t,r)(x,t)). (3.3) 

The Lagrangian approach to the CH equation is ideally suited to the multisym- 
plectic variational theory, and we begin by specifying our fiber bundle ttxy ■ Y — > 



See [FF 1981 and [CH 1993 
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X. Let X = S 1 x R, and Y = S 1 x R x R. We coordinatize X by (a; 1 , a ) (or 
(x, t)) and y by (x 1 , x°, y) (or (x, t,y)). A smooth section e C°°(F) represents a 
physical field and is expressed in local coordinates by (x, t, T)(x, i)), where 77 is the 
Lagrangian flow solving ( [3.3] ). The material or Lagrangian velocity (d/dt)rj(x,t) is 
an element of T^r x WY = Tr x j >y \Y, where y = r)(x, t). 

Using ( |3.3| ) together with u y — rjtx/Vxi the Lagrangian representation for the 
action may be expressed as 



x 



jfel 2 + ri x 1 rfi x )dxdt 



(3.4) 



The second jet bundle J 2 Y is a nine-dimensional manifold and 2-holonomic sections 
of J 2 Y — > X have local coordinates 



j 2 (4>) = (x,t 1 r)(x,t),r] x (x,t) 1 r] t (x 1 t),r] xx (x,t),r] xt (x,t),r]t x {x,t),r]tt{x,t)) , 

where for smooth sections r) x t{x, t) = ri tx (x, t). The Lagrangian density C : J 2 Y — > 
A 2 (X) is expressed as 

C(x 1 ,x ,y,y 1 ,y ,y 11 ,y 10 ,y Q1 ,y 00 ) = L^ 1 ,x° ,y,y 1 ,y ,y 11 ,y 10 ,y Q1 ,yo )dx 1 A dx° . 

For the Camassa-Holm equation the Lagrangian density evaluated along the second 
jet of a section <fi is given by 



dx A dt . 



(3.5) 



As our Lagrangian ( |3.5| ) depends only on yi,y$, and 7/01 > the Euler-Lagrange 
equation (2.14) simply becomes 



d ( dL \ d ( dL \ d 2 ( dL 



dx \ di] x ) dt \ di] t ) dtdx \ drj tx 



0, 



so that we have the Lagrangian version of the CH equation (3.1) given by 

2 \ 



^ -rt) -ta) t+ ^ =0. 



(3.6) 



(3.7) 



By differentiat ing u = (d/d^rjor]^ 1 three times, one may verify that (3.7) is indeed 
equivalent to (3.1). 

Now, using ( 2.11 ) we have that the Cartan form G>£ is given by 



dL , , dL n dL 
£ = T—d V A dt - — - D x — 
dr\ x \ dr\ t \ dr\ tx 

T dL dL dL 

L - — % - 77— »?t - q — ?7te + -Da 
orya; dryt dr\ tx 



dL , 

(277 A dx + — dr] t A dt 



dL 
dr\ tx 



di] t . 
rjt) dx Adt , 



(3.8) 



dL , , , , , , / dL n f dL 

Oc = -tt— (dr) A dt - r) x dx A dt) + — D x - — 

dr\ x \ dr\ t \ dr\ tx 

dL 

+ — (dr/t A dt — r/txdx A dt) + Ldx A dt 

drjtx 

if written in terms of the system of contact forms. 



(— drj A dx — rjtdx A dt) 
(3.9) 



IS 
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3.2. The multisymplectic form formula for the CH equation. Marsden 
Patrick, and Shkoller in their [ MPS 199S| paper have demonstrated how the multi 



symplectic form formula for first-order field theories when applied to nonlinear wave 



equations generalizes the notion of symplecticity given by Bridges in [B 1997 . Us- 
ing the example of the CH equation, we present below a simple interpretation of 
the multisymplectic form formula for the second-order field theories. We show that 
the MFF formula is an intrinsic generalization of the conservation law analogous 



to the one in Appendix D of [B 1997| 



Bridges has introduced the notion of a Hamiltonian system on a multi-symplectic 
structure. A multi-symplectic structure (Ai,u x , . . . , cu n ,oj°) consists of a manifold 
Ai, the phase space, and a family of pre-symplectic forms. The phase space Ai is a 
manifold modeled on R n+1 . A Hamiltonian system on a multi-symplectic structure 
is then represented symbolically by (Ai, u> , ■ ■ ■ , u> n , to , H) with governing equation 

87 87 87 

^(-^,v) + ■ ■ . + ^( — ,v)+u°( — ,v) = (VH(Z),v) (3.10) 

for all vector fields v on Ai where (•, •} is an inner product on TAi and Z(x 1 , . . . x n , t) 
is a curve in Ai. Bridges has shown that this formulation is natural for studying 
wave propagation in open systems. Bridges, in particular, has obtained the follow- 
ing conservation law in the case of the wave equation: 

— u°(Z t ,Z x ) + — u 1 (Z t ,Z x ) = Q . (3.11) 

This law generalizes the notion of symplecticity of classical mechanics. 

Let us make an appropriate choice of the phase space Ai for the Camassa- 
Holm equ ation. Our choice is entir ely governed by the coefficients in the Cartan 
form (|3.8| ). Since the Lagrangian (3J)) does not explicitly depend on time and 



space variables, that is the system is autonomous, we identify sections </) of Y with 
mappings n(x, t) from R 2 into R , an d similarly, sections of J 3 Y with mappings from 
R 2 into R 15 . The Cartan form ( |3.§| ) suggests to introduce the following momenta: 

(3.12) 



x _ dL t _ PL _ jj I dL \ 
" 9>7x " dm x \9»?tx J 
^tx dL pXx ptt pXt q 



k 1 drjta; 

Since Qc is horizontal over J X Y, the covariant configuration bundle is really J l Y — > 
X, and one should think of (77, r] x , rj t ) as field variables with each field variable hav- 
ing conjugate multi- momenta. For example, p x ,p t function as conjugate spatial 
and temporal momenta for the field component 77. Then the transformation 

{V,Vx,Vt,Vxx,Vxt,Vtx,Vtt, ■■■)' — > {T],i] x ,i] t ,p x ,p t ,p tx ) 

defines a mapping from the space of vertical sections of J 3 Y — > X into the phase 
space Ai = R 6 modeled over X = R 2 . We denote this transformation by FL. Let 
us now state the result that connects our paper to Bridges' theory: 

Proposition 3.1. The multisymplectic form formula (MFF) yields a multi-symplectic 
structure (Ai, u> 1 such that the MFF formula becomes an intrinsic generaliza- 
tion of the following conservation law: for any V,W in T that are -n xy -vertical 

8 8 

—u\T¥L-j 3 (V),T¥L-f(W)) + —u (T¥L-f(V),T¥L-f{W)) = . (3.13) 
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Moreover, the CH equation in both Eulerian form \3.\ ) and Lagrangian form ftS.']) is 
equivalent to the Hamiltonian system of equations on the multi-symplectic structure 
with the Hamiltonian defined by 

H = L-p x rj x -p t r lt -p tx rj tx . (3.14) 



Proof. Consider two 7rxy-vertical vectors V and W in T . Then TFL ■ j 3 (V) and 
TFL ■ j 3 (W) are vertical-over-X vector fields onlxM — > X whose components 
we shall denote via 

(v v ,v^,v T,t ,v pX ,v p \v pt *) , 

or just numerate by (V 1 , V 2 , V 3 , V 4 , V 5 , V 6 ). Thinking of the components of the 
transformation FL as functions on J 3 Y, we immediately see that 

VP" = dp x ■ j 3 {V) = j 3 {V)[p x ] 

VP' = dp 1 ■ j 3 (V) EE j 3 (V)[p*] 

yp tx = dp tx . j3 (yj - j3 (y) [pte] 

Using expressions ( 3.1 2| ) we express flc as: 

= dp x AdrjAdt — dp 1 A dn A dx + dp tx A drj t A dt 
- r] x dp x A dx A dt - rjtdp 1 A dx A dt — n tx dp tx A dx A dt . 



(3.15) 



Next, combining with (3.15) we obtain that 

f(W) J j 3 (V) J Qc — {V^W 71 - W^V* + V ptx W Vt 
- {V pt W^ - W pt V n }dx , 

so that 



W p V^dt 



3 



dU x 



o<b x L y[j 3 (w)jj 3 (v)jn c ]= (3.i6) 

f {V^W 1 - W^V 1 + V 6 W 3 - W 6 V 3 )dt - (y 5 ^ 1 - W^V^dx . 

JdUx 

The integral on the right-hand-side of the above equation leads us to introduce two 



Bi 



skew- symmetric matrices B\, Bq on 


U 6 : 























1 





" 
















1 






























































1 


B = 






















-1 

























































-1 

























-1 
































given by u> u {u, v) 



(B v u,v) 



To each matrix B u , we associate the 2-form lo v on ] 

v T B,yU, where u 7 v G R 6 . Wi th th e definition of uj u and the use of (3.16), the 
multisymplectic form formula ( [2.27 ) becomes, for Ux C X, 

w 1 (TFX • j 3 (V),TFL ■ j 3 (W))dt - u>°(TFL ■ j 3 {V),TFL ■ j 3 {W)) = . 

dUx 

Hence by the Stoke's theorem, 
d 



dUx 



dx "\TFL ■ j 3 (V), TFL ■ j 3 (W)) + —u>°(TFL ■ j 3 (V), TFL ■ j 3 (W)) 



dx A dt = 



Since Ux is arbitrary, we obtain the desired conservation law ( |3.13 ). 
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In the special case, when the components V n = r\ x and W n = rjt, one may verify 
that 

TFL-j 3 (V) = (rj, Vx , Vt ,p x ,p t ,p tx ), x , 
T¥L ■ j 3 (W) — (rj,r] x ,i] t ,p x ,p t ,p tx ) >t , 

so that, letting Z denote an element (77, r\ x , rjt, p x , p l , p tx ) G Ad, the formula (3.13) 
takes the special form 

—cj 1 (Z t ,Z x ) + —u°(Z t ,Z x ) = 0, 



which is the complete analogue of Bridges' conservation law ( 3.11 ) for the wave 
equation. 

Next, since the inner product (■, •) is independent of Z e Ad, the Hamiltonian 
system of equations ( 3.10 ) on the multi-symplectic structure (AdjUi 1 ,ui°) may be 
written as 

ZtJuj 



Z x J ui 1 
which results in 



d_ 

dx^ 



VH , or 


B\Z X 


9 t 


dH 


+ d/ = 


dr] 


d tx 


dH 


p tx = 

OX 


drjt 


d 


dH 


dx 


QpX 


d 


dH 


m q = 


dp* 


d 


dH 


a~ T l t = 
ox 


Qptx 



B Z t = VH , 



With the choice of the Hamiltonian ( [3. 14 ) the last four equations yield identities, 
and the first equation becomes 



d_ 

dx^ 



d_ 

at 



p* = . 



Using the Legendre transformation expressions (3.12) for p x ,p ', the latter equation 
recovers the Euler-Lagrange equation (3.6), hence (3.7). In other words, the Euler- 
Lagrange equations on J 3 Y are equivalent to Hamilton's equations on the multi- 
symplectic structure (Ad , lo 1 , lu° , H) . □ 

4. Discrete second-order multisymplectic field theory 

4.1. A general construction. We shall now generalize the Veselov-type dis- 
cretization of first-order field theory given in [MPS 1998] to second-order field 
theories, using the Camassa-Holm equation as our example. We discretize X by 
Z x Z = and the fiber bundle Y by X x R. Elements of Y over the base point 

(i,j) are written as j/y and the projection ttxy acts on Y by TTxyiyij) = (hj)- The 
fiber over G X is denoted Y^. 

For the general case of a second-order Lagrangian one must define the discrete 
second jet bundle of Y , and this discretization depends on how one chooses to ap- 
proximate the partial derivatives of the field. For example, using central differencing 
and a fixed timestep k and spacestep h, we have that 
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_ Vi+ij - Vi-ij _ Vij+i ~ Vij-i _ Vi-ij - 2 yiJ +y i+lj 

Vx ~ 2h ' Vt ~ 2k ' Vxx ~ 

„ ^ Vi+ij+i - Vt+ij-i + Vi-ij-i - Vi-ij+i n ^ Vij-i - tyjj+yij+i . . 
Vtx m , Vtt p , (4.1) 

where = r](xi,tj) and {(xi,tj)} form a uniform grid in continuous spacetime. 




Figure 4.1. Equivalent computational grid in the physical domain. 

We observe that an 9-tuple 

{Ui-ij-i, Vi-ij, Vi-ij+i, Vij-i, Vij, Vij+i, Vi+ij-i, Vi+ij, Vi+ij+i) 
is sufficient to approximate j 2 (j)(P), where P is in the center of the cell 
ffl« = ( (i - l,j - 1), (i - l,j), (i 1), 

(i + l.j-l), (i + l.j), (i + l.j + l)). 

Let A ffl denote the set of cells, i.e. A ffl = {ffl^ | G X}. Components of a 

cell are called vertices, and are numbered from first to ninth. A point E X 

is touched by a cell if it is a vertex of that cell. If U C V, then G X is an 
interior point of U if U contains all cells touching The interior int U of U 

is the set of all interior points of U. The closure cl U of U is the union of all cells 
touching interior points of U. A boundary point of U is a point in U and c\U 
which is not an interior point. The boundary of U is the set of boundary points, 
so that dU=(UD cl U) \ int U. 

A section of the configuration bundle Y — > A is a map <j> : U C X — ► F such 
that itxy °4> = id[/ . We are now ready to define the discrete multisymplcctic phase 
space. 
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Definition 4.1. The discrete second jet bundle ofY is given by 

J 2 Y = {{Vi-lj-li Vi-lji Vi-lj+li Vij-1) Viji Vij+lt Vi+lj-1) Ui+lji Vi+lj+l) I 

G X, yi-ij-i, ... , Vi+ij+i G R} 

EE X K 9 . 

The fiber over (i,j) 6 A is denoted J 2 Yy. We define the second jet extension 
of a section <f> to be the map j 2 4> : X — ► J 2 Y given by 

j 2 4>(i,j) = (ffly, 0(t-l,i-l), 0(t-l,j), + 

+ 0(t + l,j), + + ). 

Given a vector held u on Y the second jet extension of v is the vector held j 2 v 
on J 2 Y" dehned by 

i 2 u(yi-ii-i, ... , Vi+ij+i) = (v(yi-ij-i), 

v(Vij-i), v(yij), v{y lJ+ i), 
v(y i+lj -i), v(y i+ ij), v(yi+i j+ i)). 

Of course, this may easily be generalized to more accurate differencing schemes 
that require more than nine grid points to dehne second partial derivatives. 

4.2. A multisymplectic-momentum algorithm for the CH equation. Re- 
stricting our attention to the CH equation and noting that its Lagrangian depends 
only on rj x ,r]t, r) tx , we may significantly simply our discretization of the second jet 
bundle J 2 Y; this will substantially reduce our calculations and simplify the expo- 
sition. 

To approximate j 2 tfi(P) we choose the forward difference evaluations ofr] Xl i] t , r/ tx : 

_ Vt+ij ~ Vij n _ yij+i - Vij 

'IX ~ , i '/t 



Vtx 



h ' k 

Vi+ij+i - Vi+ij - Vij+i + Vij 



hk 

For this particular choice, our cell reduces to a rectangle. A rectangle □ of X 
is an ordered 4-tuple of the form 

Uij = + + + (*,j + l)). 

For each rectangle, D 1 , D 2 , D 3 , D 4 stand for the first, second, third, and fourth 
vertices respectively. If is the first vertex, we shall denote the rectangle by 
□y. The set of all rectangles in X is denoted by A D . The set-theoretical definitions 
of Subsection 4.1 apply here. For example, a point P = (i,j) G X is touched (see 
Fig. 4.2) by four rectangles Ely, Di-ij, Dy-i, etc. 

Again as (3.5) does not depend on rj xx ,rjtt, we may restrict ourselves to a sub- 
bundle B of the continuous J 2 Y defined via B = {s € J 2 y | s MM = for /i = 1,0}. 
Then the discrete analogue B (see Fig. 4.3) of B is identified with 

& = {(?/«•> 2/t+iiiy»+ij+i,2/ii+i) I G X, Vij, Vi+ij, Vi+ij+i, Vij+i eR} = X n xl 

For a section <j> : U C A — > F, we define the second jet extension of </> to B 

to be the map j 2 <f) : U C A — ► B via 

j 2 0(z,j) = (□«, ^(D 1 ), 0(D 2 ), 0(D 3 ), ^(D 4 )). 
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i+1 



i-1 



FIGURE 4.2. The rectangles which touch (i, j). 




Figure 4.3. Interpretation of an element of J 2 Y when X is discrete. 

Given a vector field v on Y we extend it to a vector field j 2 v on B by 

j 2 v{yij,yi + ij,y i+ ij + i,yij + i) = {v(y ij ),v(y i+lj ),v(y i+lj+1 ),v(y ij+1 )) . 

A discrete Lagrangian on B is then a function L : B — ► R of five variables 
□ij, y2 5 J/3, J/4, where the y- variables are labeled in the order they appear in a 
4-tuple. Let U be a regular subset of X, i.e. U is exactly the union of its interior 
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and boundary. Let C\j denote the set of sections of Y on U, so Cjj is the manifold 



Definition 4.2. The discrete action is a real valued function on Cu defined by 
the rule 

S(<t>)= ]T Lofcfi(i,j). (4.2) 

□ CC/;(i.j) = D 1 

Given a section <f> on U acting as 4>{i,j) = yij, one can define an element V G 
T^Cu to be a map V : U — ► TY acting as V(i, j) = (<j)(i,j), Uy), where «y is 
thought as a vector emanating from t/y = <j)(i,j). Given an element V S T^Cjj one 
can always extend it to a vector held v on Y. On the other hand, given a vector 
field v on y , V = v o is an element of T^Cy. Thus, it is sufficient to work with 
vector fields v on Y alone. 

If v is a vector field on Y, consider its restriction v \ Yi . to the fiber Yy. Let 
F% : Yij — ► Yij be the flow of v \ Y ... Then by definition of the flow, 



v(<j>{i,j)) = ^ 



Therefore there is the f-parameter family of sections on U defined by (f) X = F% o <f> 
such that cjp = (j) and ^| A _ <p x — v o </) — V . Thus, the variational principle is 
to seek those sections 6 for which 



d 

for all vector fields v on Y. 



5(2^00=0 (4.3) 



4.3. The discrete Euler-Lagrange equations. With our choice of B, the dis- 
crete Lagrangian for the Camassa-Holm equation is 

T , 1 (vi - v\ (va - Vi) 2 , h (y 3 - 2/2 - 2/4 + yi) 2 

L[yi, 2/2, 2/3, 2/4J = r 



2 \ h k 2 2/2 — 2/1 /i 2 fc 2 

(4-4) 

The variational principle yields the discrete Euler-Lagrange field equations 

(DEL equations) as follows. Choose an arbitrary point € U. Henceforth, with 



a slight abuse of notation, we shall write yy for (f>(i,j). The action (4.2), written 
with its summands containing yy explicitly, is (see Figure 3.2) 

S = h L(yij,y i+ ij,y i+ ij + i,yij + i) + Lfa-ij, yy, yy+i, Vi-ij+x) 

+ L(yi-lj-l,yij-l,Vij,Vi-lj) + L(yij-l,Vi+lj-l,Vi+lj,Vij) h — • 

Differentiating with respect to yy yields the DEL equations: 

dL. , 3L. 

5— (yy , 2/i+ij , 2/i+ij+i , 2/ij +1 ) + tt- {Vi-ij , yij , 2/y+i , 2/i-ii+i) 

dL dL 

+ -5— W-lj -UVij -l, Vij, Vi-lj) + ^—{yij-l,yi+lj-l,yi+lj,Vij) = 

0*2/3 oyi 

for all (i, j) G int {/. Equivalently, these equations may be written as 

^(D 1 ), 0(D 2 ), 0(D 3 ), 0(D 4 )) = (4.5) 



x ^ dL 

^ dyi 



i;a-(i,j)=a' 
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for all S int U. Computing and evaluating along rectangles touching an 
interior point and substituting these expressions into ([Of), we obtain the 

discrete Eulcr-Lagrange equations for the CH equation: 

(Afct/j+ij - A k y tj ) 2 (A k y i:j - Afct/i_y) 2 (Afcj/y) 2 (A fe yi_y) 2 



2hk 2 {A hyiJ ) 2 2hk 2 (A h y l - lj ) 2 2hk 2 2hk 2 

{AkVi+ij - AkVij) (Afej/ij - AkVi-ij) (AkVi+ij-i - Akyij -1) 



hk 2 (A h yij) 



hk 2 (A h y i _ 1 j) 



hk 2 {A h y ij _ 1 



(4.6) 



+ 



(AkVij-i - AkVi-xj-x) (A h y ij )(A k y ij ) (A h y lj ^ 1 )(A k y lj ^ 1 ) 



hk 2 (A h y 



i-ij-i) 



hk 2 



hk 2 



= 0. 



where 



AuVij = Vij+i - Vij, AkVi-ij = Vi-ij+i ~ Vi-ij , 
AhVij = Vi+ij - Vij, AkVi+ij = Vi+ij+x - Vi+ij ■ 



To see that (4.6) are indeed approximating the continuous Euler-Lagrange equa- 
tion (3/7), notice that the first two terms combine to approximate 




A k y. 



A k yij 



hk 



A k yjj - Akyi-ij 
hk 



A h yjj 
h 

AhVi-ij 



As to the third and fourth terms of (4.6), 



1 1 

2 h 



A k yij 



A k y 



i-lj 



Next, the fifth, sixth, seventh, and eighth terms combine as 



rjtx 

Tlx 



1 

hk 



A k y 



i+lj 



A, 



kVij 



'A 



hVi 3 



A k y 



hk 

Afc^j-i / A h y ij _ 1 



A k yjj - AkVi-ij I Aw\H- 



hk 



h 



AkVij-i - A fc ^_ij. 
hk 



' A h y. 



i-lj-l 



hk I h 

Finally, the last two terms of ( |4.6| ) approximate 

1 / A h y %j A k y i:j A/,y^-_i A fe y y _i 



The numerical scheme ( 4.6) proceeds as follows: Suppose that 

AhVij, AhVi-ij, AhVi-ij-i, AhVij-i, AkUij-i, AkVi-ij-i, AkVi+tj-i 



as known from the two previous time steps; then (4.6) may be written as 

F(A k yij, A k y i+lj , A k y i _ lj )=0. 

These are implicit equations which must be solved for j/y+i, 1 < i < N, where N 
is the size of the spatial grid. 
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4.4. The discrete Cartan form. We consider arbitrary variations which are in 
no way constrained on the boundary dll . For each G dU there is at least one 
rectangle in U touching since (i,j) £ cl U and U is regular. On the other hand, 
not all four rectangles touching are in U since int U. Therefore, each 

£ dU occurs as the I th vertex for either one, two, or three of the I e 1, 2,3,4 
and the corresponding I th boundary expressions are given by 

oyi 

■^—(yi-lj,yij, Vij+l, Vi-lj+l)V(i,j) > 

oyi 

■w—(yi-ij-i,Vij-i,yij,yi-ij)V(i,j), 
<jyz 

dL 

— (y ll - 1 ,y l+lj - U y. l+lj ,y lj )V(i,j) , (4.7) 

where y%j = 4>(i,j). The sum of all such terms is the contribution to dS from the 
boundary dU. We thus define the four 1-forms on6C J 2 Y by 

®l(yij,yi+lj, Vi+lj+l, Vij+l) ■ (%« , , V yt+lj+1 , Uyy+i) 

dL 

-{yij,yi+ij,yi+ij+i,yij+i) ■ (v yij , o, o, o), 



dyi 



®L(Vij,Vi+lj, Vi+lj+l, Vij+l) • (%<j i u v i , ' ' ». 

<9L 



9j/2 



{yij,Vi+ij, yi+ij+i, yn+i) ■ (o, Uj/i+ii) °) °): 



®L(Vij,Vi+lj> Vi+lj+l, Vij+l) ■ K« ) > i 

<9L 



= g^(Vij,Vi+lj, Vi+lj+l, Vij+l) ■ (0, 0, + 0), 

@Z,(2/ijjyi+lj!yi+li+liyy+l) ■ 7 > Vyi+ij+i , %y+i) 

= g-^{yij,Vi+ij, Vi+ij+i, Vij+i) ■ (o, o, o, 

We regard the 4-tuple (6^, 0|, 0|, 0|_) as being the discrete analogue of the mul- 
tisymplectic form 0£. Give n a vector field v on Y such that V — v o (f> 7 the first 
expression from the list (4/7) becomes [(j 2 (f>)* (j 2 v J 0^)]((i, j)), the others written 
similarly. With this notation, dS may be expressed as 



/ 



dSffl ■ V = 

(ij')eint U 



(i,j)eau 



\ 



J2 [(ftnfvje^in 1 ) 



I nC!7;i; 



J 



J2 [(j 2 <t>nfvje i L )}(tf) 



(4.8) 



/ 



\(i,j)='d ! 

4.5. The discrete multisymplectic form formula. For a rectangle □ in X, 
define the projection n\j ■ Cjj — ► B by 

tt d (0) = (□, ^n 1 ), 0(n 2 ), (f)(n 3 ), ^(n 4 )) . 
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Calculating the form n^Q L on Cu gives 



OL 



(n* a e l L )(t) ■ v = g^-Wa 1 ), 0(D 2 ), ^(D 3 ), <MD 4 )) V(D l ) 



This immediately implies that the variation (4.8) can be written as 



dStf) ■ v = 

eint U 



( 



i Dec/;!; 



+ E 

(i,j)edU 



Pelf* 

\(i,J)=p« 



(4.9) 



Define the 1-forms a\ and «2 on the s pace of sections Cjj to be the first and the 
second terms on the right hand side of (4.9), respectively. 

As in the previous section we would like to derive the discrete analogue of sym- 
plecticity of the flow in mechanics. Let 4> x be a curve of solutions of ( |4.5| ) that 
passes through at zero with V = ^x\ x _ ' 
and each A, the following holds: 



Then for each interior point 



E 



dL 



h x {U x ), A (D 2 ), ^ A (D 3 ), <^(D 4 )) = . 



i-,a-(i,j)=a l 

Differentiating these equations with respect to A at A = 0, we obtain the 

Definition 4.3. If <fi is a solution of the discrete Euler- Lagrange equations 

then a first-variation equation solution at (p is a vector V € T^Cu such that for 

each (i,j) G int U, 



E E 



d 2 L 
dykdy. 



(^n 1 ), </>(n 2 ), 0(n 3 ), ^(D 4 )) v(n k ) = o . 



(4.10) 



i;D;(ij)=D« fc=l 

By d efin ition of the forms a.\ and ct2, dS = cui+o^. Since d 2 S = 0, dai+dct2 = 0. 
Using ( fi~9| ) and denoting the vertices of □ by yi, J/2) J/3> J/4j we have that <d l L = 
j^dyi, which implies that for all I — 1,2,3,4 



E 

fc=i 



d 2 L 
dykdy I 



dy k A dyi 



Therefore, 



n* a n l L (cb)(V,W) = ^(ttd (<£))(?>□ • U,T rr D • W) 

= ^WD 1 ), • • • , 0(n 4 )) • ((up 1 ), • • • , u(n 4 )), (^(n 1 ), • • ■ , win'))) 
= E WDl) ' 0(Q2) ' ^ (d3) ' ^K^Wn 1 ) - ^(□^(□ fc )} • 

(4.11) 
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Substitution of (4.11) into the exterior derivative of the right hand side of (4.9) 
yields 



(i,j)eint U 



d 2 L 



,— \ dy k dyi 

UCV;l; fc=l 



E E 



(^(D 1 ), • • • , ^.(D 4 )) (v(D k )W(D l ) - V(n l )W(n k )) 



and 



□ cc/;(: k=l " K yl 



(i,j)edu 



(^D 1 ), • ■ • , 0(D 4 )) (V(D*)W(n<) - V(n')W(n fc )} 



Vui_u ;t ; 
(i,j)=P< 



When specialized to two first-variation solutions V and at 0, da>i (<f>) ( V, W) 
vanishes, because for each interior point all four rectangles touching it are 

contained in U, and V(D l ) = V(i,j) and W(D l ) = W(i,j). Therefore, dai = 
and the equation d 2 S — becomes da 2 = 0, which in turn is equivalent to 

( \ 

E 

(i,j)edu 



. DCC/;!; 







(4.12) 



/ 



for all vector fields v,w on Y. This is the discrete analogue of the multisymplectic 
form formula for the continuous spacetime. 

We observe that dL = <d\ + 6| + 6f + 9f , which shows that 

nj, + o| + n| + n| = o , 

which in turn implies that only three of the 2-forms Sl l L , I = 1,2,3,4 are in fact 
independent. In addition, this implies that for a given and fixed rectangle □, 

4 

o = 5>&n' £ (4>)(v,w) 



i=i 

4 4 

EE 



d 2 L 

dykdy, 



^(n 1 ), • • • , <hd 4 )) (v(n k )w(o l ) - v(n l )w(n k j\ 



i=i k=i 

for all sections <j> and all vectors V, W. 

4.6. The discrete Noether theorem. We would like to derive the discrete ver- 
sion of the Noether theorem for second order field theories. This is not the most 
general form possible as we are working with a particular example. However it is 
such as to facilitate the derivation of any other case without significant effort. 

Suppose that a Lie group G with a Lie algebra q acts on Y by vertical symmetries 
such that the Lagrangian L is invariant under the action. Vertical action simply 
means that the base elements from X are not altered under the action, hence the 
action restricts to each fiber of Y. Let $ : G x Y — ► Y denote the action of G on 
Y. For every g G G let $ g : Y — ► Y be given by y^ i — > $(g,yij). We also use 
the notation g • y — & g (y) for the action. Then there is an induced action of G on 
B defined in a natural way: 

9 ■ (2/1,2/2,1/3,2/4) = ($(#, l/i), $(t7, 2/2), $(g, 1/3), ®(g, 2/4)) • 
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Recall that the infinitesimal generator of an action (of a Lie group G on a manifold 
M) corresponding to a Lie algebra element £ s g is the vector field £m on M ob- 
tained by differentiating the action with respect to g at the identity in the direction 
£. By the chain rule, 



[exp(iO 



where exp is the Lie algebra exponential map. 
Using this formula, we immediately see that 

63(2/1,3/2,2/3,2/4) = (M2/i),M2/2),M2/3),6'(2/4)) 

The invariance of the Lagrangian under the action implies that 

£b -I dL = for all ^eg, 

which, for a given □, is equivalent to 

2=1 



^^—(2/1,2/2,2/3,^4)^(2/0 = ( 4 - 13 ) 



for all £ G g and all (2/1,2/2, 2/3, 2/4) e £>■ F° r eacn l> l et us denote by n a i : B — ► Y 
the projection onto the I component. Using this projection the four components 
of the infinitesimal generator £g are expressed as 

4 4 d 

;=1 i=l y 



Hence, the equation ( [4. 13| ) becomes 

4 

^(gje^O for all £e . (4.14) 
(=i 

We observe that for each I, 

dL 

02// 

is a function on £> which we denote by J l {£). Notice that J z (£) =(gJ 6^ is the 
discrete multisymplectic analogue of £m Jwi = dJ(£) in classical mechanics so that 
£m is the global Hamiltonian vector field of J(£). Many symmetry groups act by 
special canonical transformations, i.e. £^ m 6l = 0, in which case J(£) = £m -I 
In such case, J(£) is uniquely defined. 

Since £g is linear in £, so are the functions J ( (£), and we can replace the Lie 
group action by a Lie algebra action £ 1 — > £g. Finally we are ready to define the 
momentum maps. 

Definition 4.4. There are four g* -valued momentum mappings J , Z = 1,2,3,4 
on B defined by 

(J i (yi,2/2,2/3,2/4),£) = ^(0(2/1,2/2,2/3,2/4) (4.15) 



/or all t; £ g and (yi, j/2, 2/3, 2/4) G where (•, •) is i/ie duality pairing. 
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The equation ( 4.14) implies that 

J 1 + J 2 + J 3 + J 4 = , 

so, as in the case of the Lagrangian 2-forms, only three of the four momenta are 
essentially distinct. 

The discrete version of the Noether theorem for second-order held theories now 
follows. Define the action of the Lie group G on Cu by 

g -<f>= $ 9 of i.e. (g-(f>)(i,j) = $(g,4>{i,j)); 

since the Lagrangian is G- invariant, then 

S{g ■<!>) = J^LojVMa 1 ) 

□ CU 



acu 



£ L^P 1 ),.- ,0(D 4 ))=5(0). 



acu 

Once again letting g — exp(i£) and differentiating with respect to t at t = 0, 
we obtain that (£,Cu -I dS)((f>) = for all 4> € Cu- One can readily verify that 
S.Cui'P) — £y ° 0, which is an element of T^Cjj- Thus, 

dS{4>) ■ (£y o 0) = (4.16) 

for all £ g q and (p £ Cu- Since iS is G-invariant, then G sends critical points of S 
to themselves, or in other words, the action restricts to the space of solutions of the 
Euler-Lagrange equations. Therefore, if is a solution, so is 0* = exp(i£) • cj>, where 
<t>° = <f> and j| _ $ = £y ° <j>- Substituting 0* into the discrete Euler-Lagrange 
equations and differentiating with respect to t at t = 0, we obtain that for any £ 
and <j>, £y o <j> is a first-variation equation solution. Using (18), ( 4.16 ) becomes 



( 

o = dS(4>) ■ (fr o 0) = £ £ 

(i,j)EdU \ DCt/;!; 

\(i,i)=o' 



\ 



(^(n 1 ), ••• ,0(n 4 )Ky° </>(□') 



/ 



E 

(i,j)eac/ 



E 



(frjeLX^n 1 ), 



□ CJ/;i; 

\(i,j)=n ( 



]T J'WQ 1 ), •■•^(□ 4 ))(0 



□ CJ/;I; 



a 4 )) 

\ 

/ 



for all </> from the solution space and all £. Thus, the discrete version of the Noether 
theorem is 

/ \ 



E 



J2 [(i^n'Kn 1 ) 



(i.j)edU □££/;!; 



= 



(4.17) 
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for all <j> from the solution space. 
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